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■ Abstract. We investigate the response function of LISA and consider the 

' adequacy of its commonly used approximation in the high-frequency range of 

the observational band. We concentrate on monochromatic binary systems, such 
f**l , as white dwarf binaries. We find that above a few mHz the approximation starts 

, becoming increasingly inaccurate. The transfer function introduces additional 

i amplitude and phase modulations in the measured signal that influence parameter 

estimation and, if not properly accounted for, lead to losses of signal-to-noise ratio. 
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PACS numbers: 04.30.-w, 04.80.Nn, 97.20.Rp, 97.80.-d 

QO 1- Introduction 

O ' 

' The Laser Interferometer Space Antenna (LISA) is a collaborative ES A/NASA 

enterprise to fly a space borne laser interferometer of arm-length 5 million km for 
the observation of gravitational waves (GW's) in the frequency window 1CP 4 Hz 
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0.1 Hz £Q. LISA is likely to observe the largest sample of galactic binary systems 
of stellar-mass compact objects, strongly dominated by double white dwarfs (WD's), 
but including also neutron stars and possibly solar-mass black holes Due to the 
very large sample of these (clcctro-magnetically very faint) objects, it is therefore 
interesting to investigate new astrophysical information that LISA will be able to 
provide, and therefore estimate the errors associated with parameter measurements. 
All the investigations that have been carried out so far regarding short-period stellar 
mass binary systems have implicitly assumed LISA to work in the so-called long 
wavelength approximation (which is equivalent to assuming that the instantaneous 
GW emission frequency is / <C 10~ 2 Hz), in other words that the LISA transfer 
function is constant EJ El El EI E| • However the vast majority of the detectable 
systems - current estimates suggest that LISA will be able to resolve ~ 10 4 sources - 
will be in the frequency band 1 mHz - 10 mHz where the frequency and location 
dependent transfer function might introduce significant effects; such effect has already 
been noted for the case of observations of black hole binary inspirals 9 . 

Here we address the effect of LISA's finite arm-length for parameter estimation 
in observations of monochromatic signals, and compare the results with those from 
previous analyses. 
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2. Parameter estimation 

The signal at the LISA Michelson detector output reads ^H] 

h(t) = D ab (t)h ab (t), (1) 
where h a b(t) describes the metric perturbation induced by GW's, and 

D ab (t) = \ [ifflqm® - imil(t)T k (t)] (2) 

is the detector response tensor. In the previous expression is the unit vector along 
the j-th arms of the interferometer (j = 1,2,3 for the three arms of LISA here we 
simply consider one detector with j=l and fc=2) and 
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isinc 
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is the LISA transfer function which tends to 1 in the low frequency limit defined by 
///* < 1; /* = c/(2nL) ~ 9.5 mHz is the LISA transfer frequency (L — 5 x 10 6 km) 
and the unit vector N a identifies the source location in the sky. We note the time 
dependence of Tj{t) due to the change in orientation of LISA during the year long 
observation time pp. 

Setting T = 1 in Eq.J2J and then substituting into Eq.Q returns the 
familiar expression of the signal at the detector output in the long wavelength 
approximation 

h(t) = F+(t)h+(t)+F*(t)h*(t), (4) 

where F + and F x are the so-called antenna beam patterns and h + (t) and h x (t) are 
the two independent polarisations. We henceforth refer to these two representations 
as h for the exact response function and hi for the response function in the long 
wavelength approximation. 

To establish at what frequency the two representations h and h\ become 
significantly different, we compute the overlap 

0= - (W , (5) 

v / PKW 

where (u\v) represents the usual inner product between two functions u and v 
Notice that in Eq (JSJ we do not maximise over any of the signal parameters. We 
have calculated O over the frequency range 10 -4 Hz - 10 -2 Hz for randomly selected 
sky locations and orientations of the orbital plane, which is identified by the unit 
vector L a along the direction of the orbital angular momentum; we have found that 
O w 0.9 at w 5 mHz (the peak of the LISA sensitivity is at 3 mHz) and the overlap 
drops rapidly below « 0.5 for / ~ 10~ 2 Hz, which suggests that at frequencies higher 
than a few mHz the low frequency approximation may no longer be appropriate; this 
has important repercussions for signal detection through matched filters, and we are 
currently investigating this issue in more detail. 

To investigate the statistical errors associated with parameter estimation we 
compute the Fisher information matrix T^, and its inverse (r- jfc ) _1 = (AA 3 AA fc ), 
which is related to the variance-covariance matrix through a straightforward 
normalisation of the off-diagonal elements We do this for both h and h\. A J 
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represents the vector of the unknown parameters, in our case the signal amplitude 
A, the constant frequency /o, the arbitrary initial phase (f>o and the four angles that 
define the source geometry with respect to the detector (in our representation we use 
as independent parameters cos On, 4>n, cosQl, 4>l)- We also evaluate the LISA angular 
resolution, defined as AQ = 2tt((Acos6 n ){A<P n ) - (Acos6>at AoV) 2 ) 1/2 - 

We investigate the differences between the expected errors computed using the 
low frequency approximation and the exact expression by computing their fractional 
difference 

X id = (X } -X)/X, (6) 

where X stands for either the mean-square error or the angular resolution and the 
subscript / indicates that the long wavelength approximation has been used. We have 
calculated Xfd over the LISA observational window for randomly selected positions 
and orientations in the sky. Here we show results for A A/ A and Afi. The results 
are reported in Figures ^and [5] and clearly show that the transfer function introduces 
characteristic modulations in phase and amplitude, c.f. Eq.(3) that help in resolving 
the parameters more accurately. 

Our preliminary results indicate that LISA's angular resolution differs by « 5% 
and « 30% at 5 mHz and 10 mHz, respectively, with respect to the long-wavelength 
approximation. Analogously, AA/A differs by a few percent and « 20%, respectively. 
It is clear that the results depend strongly on the signal frequency and the source 
geometry, and hi becomes indeed an increasingly poorer representation of the actual 
detector output for / £ 5 mHz. 
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Figure 1. Comparison of the measurement errors. The plot shows the 
fractional difference of AA/A, Eq.J^J, as a function of frequency between the 
low-frequency approximation and the exact expression of the signal measured at 
the detector output for three randomly selected positions and orientations of a 
source in the sky: solid line,2. 2,4.6,2. 5,4.6; dashed line,0.90,l. 2, 1.7,1.1; dot-dashed 
line, 2. 8, 4. 0,0. 89, 3. 9; angles are listed corresponding to 6n , <j>N ; 0L > <t>L ■ 
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Figure 2. Comparison of the measurement errors. The plot shows the fractional 
difference of AQ, for the same positions and orientations in the sky, c.f. Figure 1. 



3. Conclusions 

We have shown that in the frequency range 1 mHz - 10 mHz, where most of the WD 
binary systems will be detected, the effect of the finite length of LISA arms can be 
significant and needs to be properly taken into account in LISA data analysis. Our 
analysis is limited in two main respects: (i) we have explored only a limited portion of 
the whole signal parameter space, and (ii) we have not allowed for intrinsic frequency 
drifts of the signal during the observation time. The analysis to address the two former 
issues is currently in progress jllj. 
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